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PREFACE 

This  report  was  prepared  as  part  of  RAND's  continuing  effort  in  the 
VELA  Sierra  Program. 

It  is  to  be  submitted  to  the  Physical  Review  for  publication. 


SUMMARY 

/  ft  KV’t  £ 

They^frays  from  a  nuclear  explosion  in  space  Compton— scatter 
electrons  near  the  surface  of  the  device  or  in  a  surrounding  material 


shield.  The  scattered  electrons  leave  the  surface  and  are  accelerated 


back  toward  it  by  the  positively  charged  matter.  Provided  they  are 
asymmetrically  distributed,  the  accelerating  electrons  radiate  an 
electromagnetic  signal.  The  electron  motions  are  analyzed,  the  electro¬ 
magnetic  signal  is  estimated,  and  its  detectability  is  discussed. 

For  a  typical  nuclear  explosion,  the  electromagnetic  signal  is 
independent  of  the  yield  and  contains  frequencies  up  to  10  to  100  mega¬ 


cycles  per  second  and  thus  will  penetrate  the  ionosphere.  Taking  into 
account  dispersion  by  the  ambient  interplanetary  plasma  os  ur^> — - 

✓ — CTtC-uii  unn/ cr,)j-  the  peak  electric  field  strength  at  a  distance  R  kilometers 

h  _  *  -  HZ-n-** 

from  the  explosion  is  /TTO  *  R  volts/meter.  The  pulse  length  is 

-tr  *  T 

<^rrir^  R  sec.  If  only  background  cosmic  noise  limits  detectability  of  the 

(f  yc*--.  < 

signal,  the  maximum  detectable  range  is  about  10  (  km.  . 


H  v  r  •  q 
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t  0.  ;  /•  '■  ^ 
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ELECTROMAGNETIC  RADIATION  FROM 
A  NUCLEAR  EXPLOSION  IN  SPACE 

I.  INTRODUCTION  ' 

During  the  technical  discussions  at  the  Geneva  conference  on  the 

nuclear  test  ban,  it  was  suggested  that  a  nuclear  explosion  in  space  might 

radiate  a  high-frequency  radio  signal.  A  number  of  different  mechanisms 

were  considered,  but  the  magnitude  and  spectrum  of  the  signal  were  only 

conjectured  at  the  time.  More  recently  a  particular  mechanism,  due  to 

the  asymmetric  emission  of  X  rays  from  the  explosion  causing  an  asymmetric 

acceleration  of  the  ambient  charges  in  the  low-density  medium  surrounding 

the  explosion,  has  been  quantitatively  considered  by  Johnson  and  Lippmann.  1 

The  maximum  detection  range  of  the  signal  which  they  estimate  is  about  a 
4 

few  times  10  kilometers  from  an  explosion  producing  one  megaton  of  X— ray 

energy;  the  frequency  of  the  signal  is  less  than  10  megacycles  per  sec. 

In  the  present  discussion  a  quite  different  mechanism  is  considered. 

2 

This  mechanism,  based  on  the  emission  of  y  rays  from  the  explosion 
rather  than  X  rays,  produces  a  signal  of  greater  magnitude  and  with  higher 
frequencies.  An  estimate  of  this  signal  will  be  given  here  and  the  possibility 
of  observing  it  will  be  discussed. 


^M.  H.  Johnson  and  B.  A.  Lippmann,  Phys.  Rev.  ,  119,  827  (1960). 

2  — 

This  mechanism  was  discussed  by  u r.  O.  1.  Leipunski  at  the  1958 

Geneva  Conference  of  Experts. 
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II.  GENERAL  CONSIDERATIONS 

During  the  course  of  a  nuclear  explosion,  y  rays  are  produced  — 
both  directly  in  the  fission  process,  and  indirectly  by  inelastic  scattering 
of  neutrons  in  the  materials  of  the  device.  These  y  rays  while  moving 
through  the  device  or  any  ambient  matter  can  Compton— scatter  electrons. 

If  the  scattered  electroi  s  are  produced  sufficiently  close  to  the  surface  of 
the  matter,  they  can  leave  it  and  start  to  move  into  the  surrounding  space. 

As  the  electrons  leave,  the  matter  becomes  positively  charged  and 
accelerates  the  electrons  back  towards  it.  Provided  they  are  asymmetrically 
distributed,  the  accelerating  electrons  radiate.  There  are  at  least  two 
possible  ways  by  which  such  an  asymmetrical  distribution  of  electrons 
might  be  generated:  (1)  because  of  the  design  of  the  device,  the  y— ray 
flux  at  the  surface  may  be  asymmetric,  or  (2)  the  y  rays  may  be  emitted 

3 

isotropically,  but  there  may  be  an  asymmetric  external  material  shield. 

It  is  the  electromagnetic  radiation  generated  by  the  Cempt on-recoil 
electrons  from  the  device  or  from  an  external  shield  which  we  wish  to 
estimate.  Some  of  the  essential  quantities  determining  this  radiation, 
however,  depend  upon  the  design  of  the  specific  device;  we  shall  use  an 
idealized  and  simplified  model  of  the  explosion  device  and  typical  values 
for  the  relevant  parameters  —  thereby  obtaining  only  an  order  of  magnitude 

3 

Such  shields  have  been  proposed  to  reduce  the  X— ray  flux  radiated 
in  a  given  direction  from  the  explosion  and  thus  make  more  difficult  long- 
range  detection  of  X  rays.  The  shield  might  be  a  flat  plate  or  a  hemisphere 
a  few  meters  from  the  device  so  oriented  as  to  shadow  the  explosion  from 
possible  X— ray  detectors. 


RM— 2840 
3 


estimate  of  the  intensity  and  detectability  of  the  radiation  generated  by  the 
Compton— recoil  mechanism. 

III.  SOURCE  OF  THE  SIGNAL 

During  the  course  of  the  nuclear  explosion  the  y  rays  will  be  assumed 
to  be  produced  at  approximately  an  exponentially  increasing  rate.  If  t  =  0 

is  the  time  at  which  the  explosion  starts,  then  the  number  of  y  rays  produced 

#  t  8  _ 

up  to  time  t  is  given  by  e  ,  where,  typically,  or  is  about  10  sec  .  As 

the  explosion  proceeds,  the  generation  of  y  rays  continues  to  increase 

exponentially  until  a  substantial  fraction  of  the  full  yield  has  been  produced, 

at  which  time  the  rate  of  generation  of  y  rays  reaches  a  maximum.  The 

generation  from  then  on  proceeds  at  a  considerably  slower  rate  than  the 

build-up.  For  simplicity,  we  shall  assume  the  generation  increases 

exponentially  at  the  rate  oreat  until  time  T  and  then  decreases  to  zero 

* 

exponentially  at  the  rate  cr/3e ~  (a  >>  0,  u  £  1).  About 

three  hundredths  per  cent  of  the  total  yield  is  assumed  to  be  produced 

as  y  rays.  These  y  rays  have  an  average  energy  of  about  1  Mev,  and 

21  19 

there  are  7.5  x  10  y  rays  produced  per  kiloton  (4  x  10  ergs)  of  yield  — 

thus, 

aT  21 

(1  +  ff)e  =  7.5  x  10  Y  ,  (1) 

where  Y  is  the  total  yield  of  the  explosion  expressed  in  kilotons. 

The  y  rays  Compton— scatter  electrons  within  the  device  or  in  a 
surrounding  shield,  if  one  is  present.  The  resultant  recoil  electrons  will 
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be  distributed  both  in  energy  and  direction,  but  will  move  predominantly 

in  a  forward  direction  with  a  mean  energy  of  about  1/4  Mev  —  corresponding 

to  those  electrons  produced  at  a  depth  of  one— half  the  electron  range  R@ 

2 

(R^  %  0.  1  g/cm  )  from  the  surface  of  the  device  or  the  shield. 

The  total  number  of  electrons  emitted  per  unit  area  up  to  time  t  is 


.  .  Re  1  art 
n(t)  =  2T"  6  ’ 

y  4ira 


for  t  <  T, 


R 


2X 


y  4?ra 


1  aT  f. 

“2  e  L 


1  +  cr  —  ere 


— f3(t— T) 


,  for  t  >  T,  (2) 


where  a  is  the  radius  of  the  nuclear  device  or  the  distance  from  the  source 

of  y  rays  to  the  surface  of  the  shield,  and  X^  is  the  mean  free  path  for 

2 

Compton  scattering  (X^  %  15  g/cm  ).  If  the  thickness  of  the  matter  in  which 
the  Compton— scatterer'  electrons  are  produced  is  less  than  R  ,  the  electron 
range,  then  Rg  in  liq,  (2)  should  be  replaced  by  the  thickness  of  the  matter 
and  the  mean  exiergy  of  the  electrons  is  increased.  The  limiting  mean 
energy  for  small  thicknesses  is  about  1/2  Mev. 

For  simplicity  we  shall  henceforth  consider  only  the  problem  of  a  thick 
hemispherical  shield  of  radius  a  centered  about  the  device.  This  problem 
also  provides  an  estimate  of  the  radiation  from  a  flat-plate  shield  or  from  a 
spherical  device  of  radius  a  whose  y— ray  asymmetry  is  roughly  equivalent 
to  a  distribution  wherein  y  rays  escape  predominantly  into  a  hemisphere.^ 

4 

More  generally,  it  is  equivalent,  in  first  order,  to  any  asymmetry 
in  which  the  electrons  are  distributed  in  angle  approximately  as 

Yo(0,  4>)  +  Y°(0,  <j>). 
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As  electrons  leave  the  surface  of  the  hemispherical  shield  it  becomes 

positively  charged  and  the  electrons  are  accelerated  back  towards  the 

shield.  Assuming  the  electrons  are  all  emitted  in  a  direction  normal  to 

5 

the  surface  (and  that  they  do  not  move  far  from  the  surface  ),  a  dipole 
moment  normal  to  the  surface  is  generated.  From  symmetry  the 
dipole  moment  is  uniform  over  the  hemisphere  and  has  a  strength  per  unit 
area  given  by 


Z(t)  =  —  e  L  x.(t)  r,(x  (t))  ,  (3) 

i  1  1 

where  q(x)  =  0  if  x  <  0,  rj(x)  =  1  if  x  >  0,  and  x^t)  is  the  distance  of 
til 

the  i  electron  from  the  surface.  The  sum  includes  all  electrons  which 
are  emitted  from  the  surface  in  the  unit  area  considered. 

The  acceleration  of  the  dipole  moment  is  foimd  by  summing  the 
accelerations  of  the  individual  electrons  outside  the  surface  and  including 
a  contribution  due  to  the  flux  of  electrons  in  and  out  of  the  surface,  namely, 

Z(t)  =  —  e  £  [  x.t)(x  )  +  v^  6(x.)  +  v^S(x  )]  ,  (4) 

j  1  1  U  1  V  1 

where  v^  and  v^  are  the  speeds  of  electrons  leaving  and  returning  to  the 
surface  and  6(x)  is  the  delta  function.  Strictly,  Eq.  (4)  should  be  averaged 
over  the  spectrum  of  electron  energies.  However,  for  the  present  order 
of  magnitude  estimates,  it  is  sufficient  to  assume  that  all  electrons  have 

5 

‘'From  the  solution  of  Eq.  (5),  it  is  easy  to  show  that  this  approxima¬ 
tion  is  valid. 
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the  same  energy,  namely,  their  mean  energy.  Quantitatively,  this  assumption 

is  unimportant.  However,  it  introduces  certain  artificialities  into  the  behavior 

of  the  dipole  moment  which  will  be  discussed  and  treated  below. 

If  we  neglect  magnetic  and  retardation  effects,  as  well  as  relativistic 
6  th 

effects,  the  acceleration  of  the  i  electron  is  given  by 

x.  =  —  —  E(x.,  t)  .  (5) 

i  m  i’ 

5 

Provided  x.  is  much  smaller  than  the  dimensions  of  the  shield  ,  the 

l  ’ 

electric  field,  E(x_,,  t),  is  determined  by  the  local  surface  charge  density 
only.  Thus 

E(x.,  t)  =  2ire  £  [  1  —  r|(x.  —  x.)]  +  27re  £  n(x.  ~  x  ) 

1  j  1  •>  3  •>  1 

=  4-rre  £  r|(x.  —  x.)  .  (6) 

3  J  1 

For  a  constant  electron  emission  rate  it  is  possible  to  obtain  an 
exact  solution  of  Eqs.  (5)  and  (6)  for  the  average  motion  of  the  electrons 
and  of  Eq.  (3)  for  the  (time— independent)  dipole  moment.  This  steady 
state  solution  is  presented  in  Appendix  1. 

To  determine  the  average  motion  of  the  electrons  for  a  time-dependent 
electron  emission  rate  requires  treatment  of  the  complicated  self-consistent- 
field  problem  expressed  by  Eqs.  (5)  and  (6).  Rather  than  attempt  an  exact 
solution,  we  will  estimate  the  solution  for  times  t  <  T  by  two  independent 

g 

All  these  effects  are  less  than  a  few  per  cent. 
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approximations  —  one  accurate  at  early  times  and  the  other  accurate  at 
late  times.  In  the  first  case,  the  screening  of  the  electrons  by  one 
another  will  be  disregarded.  That  is,  each  electron  will  be  assumed  to 
move  in  a  constan1  lerating  field  equal  to  the  field  at  the  emission 
surface  at  the  me  electron's  emission.  At  the  very  early  times 

during  wh  -litted  electrons  are  moving  only  outwardly,  the  no- 

screening  ap,,  .^nation  is  exact.  As  electrons  emitted  at  an  earlier 
time  begin  to  return  to  the  emitting  surface,  screening  becomes  important 
and  the  no— screening  approximation  provides  only  a  semi— quantitative 
description  of  the  electron  behavior.  In  order  to  obtain  a  more  accurate 
description  for  the  late— time  behavior  of  the  electrons,  a  second  and 
independent  approximation  will  be  made.  Namely,  the  electron  motions 
will  be  assumed  to  be  quasi-static.  That  is,  at  each  instant  the  electrons 
will  be  assumed  to  be  in  that  collective  steady  state  motion  corresponding 
to  the  instantaneous  value  of  the  emission  rate.  The  validity  of  this 
approximation  will  be  discussed  below.  Finally,  with  the  aid  of  the  results 
derived  for  t  <  T,  the  solution  for  t  >  T  will  be  obtained. 


1.  No— screening  Approximation  —  t  <  T  . 

When  screening  is  neglected,  Eqs.  (5)  and  (6)  reduce  to 


*i(t>  =  ^  f  "W 


— — — — N(t.) , 
m  i  ’ 


(7) 


where  N(t  J  is  the  number  of  electrons  per  unit  area  outside  the  surface  at 
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time  t.  ,  which  is  the  time  that  the  i  electron  is  emitted.  The  motion 
is  parabolic,  and  the  i*  electron  returns  with  the  initial  speed  Vq  at  the 
time 


t  (t.)  =  t.  + 
r  l  l 


2ire  N(t.) 


For  times  before  the  first  electron  returns, 


(8) 


N(t)  =  n(t),  t  <  tQ  +  |  ,  (9) 

and  for  times  after  the  first  electron  returns, 

N(t)  =  n(t)-/  dr  t>t  +  ±,  (10) 

r 

where  n(t)  is  given  by  Eq.  (2).  Here  tQ  is  the  time  at  which  the  first 
returning  electron  is  emitted,  and  is  found  by  differentiating  Eq.  (8)  with 
N(t)  =  n(t): 

^0  „  2  2Xv  mvo“ 

e  =  4xa  - -  .  (11) 

e  2?re 


Eqs.  (8),  (9),  and  (10)  cannot,  in  general,  be  solved  analytically  for  N(t), 
but  can  be  solved  numerically.  In  Fig.  (1)  we  show  N(t)  versus  ff(t  —  t^) 
obtained  numerically  for  early  times. 

The  qualitative  behavior  of  N(t)  for  early  times  is  simply  explained. 
Until  time  t  +  1/a  no  electrons  have  returned,  so  N(t)  just  increases 
exponentially  as  n(t).  At  t^  +  1/a  ,  there  is  a  sudden  flux  of  returning 
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electrons  which  momentarily  exceeds  the  emission  rate,  causing  N(t)  to 
decrease.  Electrons  emitted  during  this  period  will  have  a  decreased 
acceleration  and  will  return  at  a  slower  rate,  which  shows  up  as  an 
increased  N(t)  at  a  later  time.  This  oscillating  behavior  persists,  but 
the  time  between  discontinuities  decreases  as  N(t)  increases.  That  N(t) 
increases  is  expected  since  the  emission  rate  is  increasing. 

The  magnitude  and  frequency  of  the  oscillations  about  the  mean  growth 
rate  of  N(t)  depend  on  t^  and  hence  on  the  initial  velocity  v^.  Since  there 
is  actually  a  distribution  of  electron  velocities  instead  of  a  single  one,  it 
is  clear  that  in  the  actual  case  the  oscillations  will  be  substantially  averaged 
out  and  only  the  smooth  growth  of  N(t)  will  remain. 

From  the  solution  for  N(t),  Eq.  (7)  can  be  integrated  for  the  electron 
motion.  Equation  (3)  then  determines  Z(t) ,  the  dipole  moment  per  unit 
area.  Namely, 

Z(t)  =  -  e  f  dr  n(T)f"vn(t-T)-  — N(T)(t-T)2].  (12) 

t  ( T  )>  t>  T  L0  “  J 

r 

For  times  before  t  =  +  1/ur,  Z(t)  is  given  simply  by: 


mvo  c(t“V 

Z(t)  =  -^e  - 


i  c(t-V 

46 


«( t  -  tQ)  < 


(13) 


Thus,  as  is  easily  shown,  Z(t)  reaches  a  maximum  and  starts  to  decrease 
before  any  electron  has  returned.  Its  maximum  occurs  at  o'(t  —  t^)  =  in  2, 
For  later  times  Z(t)  is  shown  in  Fig.  2  where  one  expects  again  that  in 


RM— 2849 
10 


the  real  case  the  oscillations  will  not  be  present. 

For  large  t ,  we  can  obtain  an  asymptotic  expression  for  N(t) 
using  the  fact  that  the  time  between  emission,  tj  ,  and  return,  t  , 
short  compared  to  1/a  and  t.  Thus 


and 


N(t)  n(t)  -n(tj) 

mvQ«  a(  t-tQ) 

=  - 

2we 


1  -  e 


navna  e(t-4  ) 

TT*'  %(t-V 

2»re 


mv„ 


2 ire  NCt^) 


-«(t-tQ) 
e _ _ 

«2(t — tj> 


Finally,  then, 

-«(t-t  )/2 

«(t  - tj)  a;  e 


and 


mv  a  a(t~i  )/ 2 

N(t)  -  - \  e  U 

2?re 


Tfyis  asymptotic  limit  for  N(t)  is  also  shown  in  Fig.  1. 

In  this  limit  we  find  for  the  dipole  moment  per  unit  area 
2  r 


mv. 


z{t)  a  - 


i  -l/2«(t-t0) 
1+  46 


by 

is  very 


(14) 


(15) 

(16) 


12jre 


(17) 


t 
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which  has  the  very  small  acceleration 


mv„  2  -a(t-t)/2 


(18) 


The  results,  Eqs.  (16),  (17),  and  (18),  for  N(t)  and  Z (t)  are  easily  seen 

ar(tH;0)/  2 

to  hold  for  times  t  such  that  e  >>  1;  and  we  note  that  the  undesired 

oscillations  have  implicitly  been  omitted. 

2.  Quasi-static  Approximation  —  t  <  T  . 

The  no— screening  approximation,  Eqs.  (16),  (17)  and  (18),  is  not 
expected  to  be  accurate  at  late  times  when  screening  becomes  important. 
However,  at  late  times  an  alternative  approximation  becomes  valid.  Namely, 
if  the  electron  emission  rate  changes  slowly  compared  to  the  time  between 
an  electron's  emission  and  return,  that  is,  if 


E(t  -t  }  =  a(\  -  t  J  <<  1, 
n  '  r  1  r  1  ’ 


then  we  may  assume  that  the  electron  motions  at  a  given  time  are  in  that 
state  of  steady  motion  corresponding  to  the  instantaneous  emission  rate. 
From  Appendix  1  we  find  the  exact  solution  for  a  steady  state  to  be 

1/2 


N 


ss 


nmv 

=  n(tr-tl)  =  ( - T) 

ire 


(19) 


Assuming  that  a(tf  —  t^)  <<  1 ,  we  obtain  a  good  estimate  of  N(t)  and 
Z(t)  by  substituting  the  slowly  changing  rate  n(t)  from  Eq.  (2)  into  Eq.  (19). 
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Thus 


mv  a  Qf('t-t  )/2 
N  (t)  ^  ^  e  °  . 


ss 


<J~2 


or(t  -t^y.  \f2  e 
2 


-df(t— tQ)/  2 


Z(t)  ~  - 


0 


8ire 


(20) 


(21) 

(22) 


From  Eq.  (21)  it  is  deduced  that  the  quasi— static  solution,  Eq.  (22),  is 

ff(t-tQ)/2 

accurate  for  late  times  t  such  that  e  >>  1  . 

A  comparison  of  Eqs.  (17)  and  (22)  shows  that  the  no-screening  and 
the  quasi— static  approximations  agree  to  within  a  factor  of  2/  3  and  both 
predict  essentially  no  radiation  at  late  times.  This  close  quantitative 
agreement  suggests  that  in  fact  the  no-screening  approximation  is  actually 
quite  good  and  can  be  expected  to  provide  a  rather  accurate  estimate  of  the 
over— all  behavior  of  the  electron  motions. 


3.  Solution  for  t  >  T. 

In  view  of  the  behavior  of  Z(t)  for  t  <  T  ,  it  is  clear  that  for  t  >  T  , 
Z(t)  will  be  small  until  t  >>  T.  Therefore,  in  order  to  estimate  Z(t)  for 
t  >  T  ,  we  can  for  convenience  arbitrarily  assume  that  at  t  =  T  the 
electrons  are  in  a  state  of  steady  motion  (using  the  results  of  the  quasi- 
static  approximation).  In  view  of  the  accuracy  of  the  no— screening  approxi¬ 
mation,  this  approximation  can  be  used  to  determine  both  the  steady  motion 
at  t  =  T  and  the  subsequent  time— dependent  motion  for  t  >  T.  Using 
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*•  2 

Eq.  (12)  and  the  results  of  Appendix  2,  it  readily  follows  that  Z(t)  ji  for 

t  >  T.  However,  from  Eqs.  (13)  and  (17),  Z(t)  ~  o2  for  t  <  T.  Since 
2 

(|3/a)  <<  1 ,  the  radiation  at  late  times,  t  >  T  ,  is  small  compared  to 

that  at  early  times,  t  <  T.  In  addition,  it  consists  of  lower  frequencies 
which,  as  can  be  shown  from  subsequent  results,  are  more  difficult  to 
detect  at  great  distances.  Accordingly,  the  radiation  at  times  t  >  T  may 
be  neglected. 

Combining  the  results  of  sections  1  —  3,  we  conclude  that  the  dipole 
moment  will  increase  almost  exponentially  to  a  maximum  value,  decay 
somewhat  more  slowly  to  a  steady  state  value,  then  increase  again  slightly 
and  finally  decay  slowly  to  zero.  The  important  times  t,  during  which 
radiation  is  emitted,  are  in  the  neighborhood  of  cRt  —  t^)  equal  to  the  order 
of  one  or  a  few. 

For  convenience  in  subsequent  order  of  magnitude  calculations  an 
approximate  expression  for  the  dipole  moment  will  be  used.  This  expression 
was  chosen  to  give  the  qualitative  behavior  of  the  dipole  moment  as  well  as 
correct  values  for  its  maximum  and  quasi— steady  state.  Namely,  we  shall 
use 

2 

mv 

Z(t)  =  -  {  1  +  tanh  [or(t  -  t^}  ]  +8  sech  a(t  —  t^) }  ,  (23) 

where  <*t^  =  nt^  +  In  2.  For  comparison  this  function  is  shown  in  Fig.  2 
along  with  the  more  exact  results. 


I 
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IV,  THE  NEAR  RADIATED  SIGNAL 

The  electromagnetic  field  radiated  by  the  Compton  electrons  is 
composed  of  a  coherent  and  an  incoherent  part.  The  coherent  part  arises 
from  the  mean— time  behavior  of  Z(t)  given  by  Eq.  (23).  The  incoherent 
part  arises  from  the  fluctuations  in  the  mean  variation  of  Z(t).  In 
Appendix  3  an  upper  limit  on  the  magnitude  of  the  incoherent  radiation 
is  determined.  Compared  to  the  coherent  part,  even  this  upper  limit  on 
the  incoherent  part  is  quite  small  and  therefore  the  incoherent  radiation 
may  be  neglected. 

The  coherent  part  of  the  electric  field  radiated  to  a  distance  R  from 
the  explosion  is^ 

E(R,  t)  = - -  f  du  Z(a))u2  f  Yd2  re  rCt  e  x  (e  x  ep)  ,  (24) 

c  R 

where 

Z(t)  =  J  du  e1U>t  Z(u)  (25) 

defines  Z(u)  .  The  second  and  third  integrations  in  Eq.  (24)  are  carried  out 
over  the  radiating  surface,  t  ^  is  the  retarded  time  from  the  observation 
point  to  the  element  of  area  at  the  point  r  .  e  is  a  unit  vector  in  the 

1  -  r  tN  W  ft 

direction  of  the  dipole  moment  at  the  point  r  and  e^  is  a  unit  vector  parallel 
to  R  .  For  a  hemisphere  of  radius  a  with  a  radiai  dipole  moment,  we  find 

7 

The  effect  of  dispersion  by  the  ambient  medium  is  treated  in  the 
next  section. 
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E(R,  t)  = 


27ra 

o.R 


/dw 


ico(t— R/  c)  £ 


2  Jl(f  Sin9! 


Z 


(26) 


where  9  is  the  angle  between  the  axis  of  symmetry  of  the  hemisphere  and 
the  direction  to  the  observer.  The  electric  field  is  polarized  in  the  plane 
containing  these  two  directions. 

For  Z(t)  given  by  Eq.  (23)  the  Fourier  transform  of  Z(t)  is 


(27) 


From  this  result  we  conclude  that  frequencies  much  larger  than  a  are 
sharply  reduced.  Moreover,  as  can  be  seen  from  Eq.  (26),  the  finite 
size  of  the  hemisphere,  expressed  by  the  factor  J^(wa  sin  6/c)/io  ,  tends 

g  _ ^ 

to  cut  off  frequencies  above  about  1/T0  =  c/a  sin  0  .  For  a  =  10  sec  , 
and  a  =  10  meters,  the  radiated  signal  will  therefore  contain  frequencies 
only  up  to  about  10  to  100  megacycles  per  second.  From  Eqs.  (23),  (26), 
and  (27),  the  maximum  signal  strength  is  about 

2 

mv  .2  , 

E  %  0  p-  (  — )  Z  10  (R.  )  volts/meters,  (28) 

max  8eR  \  c  /  km 

2  _7 

for  1/2  mvg  =  1/4  Mev;  the  signal  duration  is  %  10  seconds. 

V.  THE  DISTANT  RADIATED  SIGNAL 

> 

To  determine  the  signal  radiated  to  great  distances  it  is  necessary 

to  take  into  account  the  fact  that  interplanetary  space  is  not  empty,  but 

2  3 

filled  with  a  low-density  ionized  plasma,  containing  from  10  to  10  electrons 


(31) 


I 
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The  integral  in  Eq.  (31)  can  be  evaluated  by  the  method  of  steepest 
descents  for  times  such  that 


“p<t_to_R/c)^  K 


(32) 


The  points  of  steepest  descent  occur  at 


V^O 


[  (t  —  t^J  —  R  / c  J 


=  +  w  . 
-  P  0 


(33) 


Noting  that  Z(— u>)  =  fZ(oo)]  ,  we  find 


•Ee|e1W 


4e‘V'-i0^0.( 


Vo}' 


(34) 


This  represents  a  rapidly  oscillating  electric  field  whose  amplitude  rises 
from  nearly  zero  at  t  =  tj  +  R/c  to  a  maximum  and  then  decays  as  t 
for  late  times.  The  maximum  amplitude,  using  the  approximation  for 

A 

Z(co)  given  by  Eq.  (27),  occurs  at  x  SQ’/w  ,  that  is,  at 


t-ti-s*  Ji  R 

0  c  2  c 
a 


(35) 


and  is  approximately 


,3/2 


3/  2  rnvn 

Emax(R)%(2,ri"  -l/-2~372  13  Jl(a  a  Sin  ®/c)  5?T 

c  R  p 


I 


(36) 
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3  5—1 

Assuming  a  %  10  cm,  5:  6  x  10  sec  (corresponding  to  an  electron 
density  of  10^/ cc),  at  %  10^  sec  *,  1/2  mv^  =  1/4  Mev,  it  is  found  that 

E  (R)  £  10  “(R.  )  ^  ^  volts/meter ,  (37) 

max  km  * 

and  the  signal  duration  is  approximately 

At  =  t  —  t'  —  R /c  J  10  ^  R.  sec.  (38) 

0  km 

7  8 

The  band  width  is  10  to  10  cycles/sec. 

VI.  DETECTABILITY  OF  THE  SIGNAL 

Since  the  frequencies  in  the  signal  extend  up  to  about  10  to  100  mega¬ 
cycles,  the  signal  can  penetrate  the  ionosphere,  and  can  be  observed  from 
the  surface  of  the  earth,  provided  background  noise  is  not  too  great.  We 
shall  estimate  the  maximum  detection  range  of  the  signal,  assuming  that 
the  only  interfering  noise  arises  from  background  cosmic  noise.  Thus  it 
will  be  assumed  that  the  receiver  can  be  effectively  shielded  against  man¬ 
made  disturbances. 

In  the  frequency  region  from  10  to  100  megacycles  per  sec  the  back¬ 
ground  cosmic  noise  is  characterized  by  an  effective  temperature  of  about 
4  o 

10  K.  That  is,  the  mean  cosmic  noise  flux  per  cycle  over  the  frequency 
range  of  the  signal  is 

— 22  T  2 

1.7  x  10  —5-  watts/m  cps 

XT 

_ 1  D 

1.7x10  „  ,  2 

- - -  watts/ m  cps. 


(39) 


I 
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If  the  detection  receiver  has  a  band  width  a  cps,  the  average  cosmic 
noise  power  received  will  be 


PN  =  = 


,8 


1. 7  x  10 


-18 


watts 


(40) 


For  a  band  width  er  it  10  cps  and  V  %  10  meters  (corresponding  to  a 

7 

frequency  of  3  x  10  cps) 


_ 11  9 

P  3M0  watts/m  . 
JN 


(41) 


The  signal  radiated  by  the  explosion  has  a  peak  intensity  of 

„2  , 

cE  ,  „ 

max  _  3  x  10  ,,  ,  2 

■  % - ;  watts/  ui  . 


4x 


(42) 


®km> 


For  distances  greater  than  10^  to  10^  km,  the  peak  signal  power  is  less 
than  the  mean  noise  flux.  For  detection,  however,  the  peak  signal  power 
need  not  exceed  the  mean  noise  power.  By  suitable  circuit  design  of 
detectors  the  mean  noise  power  can  be  biased  out  leaving  only  the  fluctua¬ 
tions  of  the  noise.  It  is  sufficient  for  detection  that  the  total  signal  energy 
exceed  a  threshold  level  determined  by  the  requirement  that  the  probability 
for  a  noise  fluctuation  to  exceed  the  threshold  is  less  than  some  acceptable 
value.  A  reasonable  threshold  for  detection  is  ten  times  the  dispersion  of 
the  total  noise  energy  during  a  time  equal  to  the  signal  duration  time. 

In  Appendix  4  it  is  shown  that  for  a  receiver  with  negligible  internal 

2 

noise,  having  an  effective  antenna  area  A  m  ,  a  frequency  band  width 
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a  eps,  and  which  integrates  the  noise  power  for  a  time  period  t  seconds, 
the  dispersion  in  the  received  energy  is 


(  -  l2" 

(£n  5n> 

1/2 

/  Or  AX.2 

\  4?rkT  .  _  ,  —22  i 

'  !7tA\ 
^  2  / 

S 

V  2r  "  , 

)  r  -0.7x10  ( 

L  J 

X. 

1/2 


T  joules.  (43) 


Clearly,  optimum  detection  results  when  a  equals  the  signal  band  width 
and  t  equals  the  signal  duration  time.  Due  to  the  dispersive  effects  of 
the  interplanetary  medium  the  signal  duration  depends  upon  distance  from 
the  explosion.  An  estimate  of  the  signal  duration  is  given  by  Eq.  (35). 

From  Eq.  (31)  the  total  received  signal  energy  is 


«  =  £  A  f  dt[E(R,  t)] 


4£¥/  k(^™)]2|zu 

cR  ui  c 


For  the  approximation  to  Z(u>)  given  by  Eq.  (27), 


E 


,  24 

2  m  vQa 


2 

e  c 


/2a  a 

\  7TC 


sin  0 


) 


(44) 


(45) 


T  (  2a  a  „\ 

l\~  BmQ) 

=  7.9, 

for 

sin  9=2.  12 

7TC 

=  12.0, 

for 

— ~  sin  0=1.  06. 
ire 

2 

./2aasin0\  ,  2aa  .  „ 

=  46.  5  ( - )  ,  for  -  sin  0  <<  1  . 

\  7TC  /  ’  7CC 


where 
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Clearly,  the  signal  energy  t  is  a  rather  insensitive  function  of  the  angle  6 , 

except  when  9  is  very  small.  Moreover,  detailed  evaluation  of  the  integral 

of  Eq.  (44)  indicates  that  nearly  all  the  signal  energy  is  contained  in  a  band 

2  7  8 

width  Ato  x;  8or /  jr,  or  u  &  4a/ ir  =  4  x  10  cps  for  a  =  10  .  Using  this  band 

2  2 
width  a  -  4  a  I  it  and  an  integration  time  r  =  (cWct)  R/c  ,  the  maximum 

detection  range  (corresponding  to  a  signal  equal  to  N  times  the  dispersion 

of  the  noise)  is  given  by 


9  2  4 

“  m  v-» 


/a\  V0  ,  „  /OrAk  \ 

(a)  -2— 1  ‘  N(-sr) 


r  4  4  8  3  2  2 

-  JL  a  m  vog 

^max  8ir  2  2.  .2  4 

L  N  w^(kT)  ce 


For  a  =  X  =  10  m,  u>  «  6  x  105  sec  1,  A  =  102  m2,  T  =  104  °  K,  N  =  10, 

’  p  ’ 

and  I  =  10  we  find 


R  =1.3x10  km . 
max 


We  note  that  the  detection  range  is  independent  of  the  explosion  yield. 

»T  at0 

This  result  is  valid  so  long  as  e  >>  e 
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VII.  DISCUSSION 

In  treating  the  propagation  of  the  explosion  signal  effects  of  the 

earth's  ionosphere  were  ignored.  To  justify  neglecting  ionospheric  effects, 

we  note  first  that  the  dispersion  of  the  signal  depends  approximately  only 

2 

on  the  total  number  of  electrons  per  cm  along  the  propagation  path.  For 
distances  greater  th^an  about  10®  to  10®  km,  the  total  number  of  ionospheric 
electrons  is  small  compared  to  the  total  number  of  interplanetary  electrons. 
Second,  we  note  that  the  cut-off  frequency  of  the  ionosphere  (~  8  megacycles 
per  sec)  is  small  compared  to  the  signal  band  width  40  megacycles 
per  sec).  Consequently,  only  a  few  per  cent  of  the  signal  energy  cannot 
penetrate  the  ionosphere. 

The  key  uncertainty  in  this  method  of  detecting  nuclear  explosion  is 
the  feasibility  of  eliminating  all  sources  of  background  noise  except  cosmic 
noise.  Removing  this  uncertainty  will  require  experimental  detector  studies. 


Fig.  2  —  Dipole  moment  per  unit  area  as  a  function  of  time 
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Appendix  1 

STEADY  STATE  MOTION 

For  a  constant  emission  rate,  the  electron  motions  approach  a  steady 
state.  In  this  steady  state  the  electric  field  distribution  is  independent  of 
time  and  is  given  by 

=  -  4jrep(x)  ,  (1.1) 


where  p(x)  is  the  electron  density  at  a  distance  x  from  the  emission 
surface.  p(x)  is  related  to  the  emission  rate  nQ  by 

p(x)  =  2  ’ 


(1.2) 


where  v(x)  is  the  velocity  of  the  outward  moving  electrons  at  x  .  (The 
inward  moving  electrons  have  an  equal-  and  opposite  velocity. )  From  the 
electron  equation  of  motion,  Eq.  (5), 


x 


=  v(x) 


dv(x) 

dx 


(1.3) 


Differentiating  Eq.  (1.  3)  with  respect  to  x  and  using  Eqs.  (1.  D  and  (1.  2), 


we  find  that 


d2v2(x)  =  16?re2  nQ 

,  2  m  v(x) 

dx 


(1.4) 


A  straightforward  solution  of  this  equation  gives  for  the  velocity  of  the 


outward  moving  electrons 
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v(x)  = 


where 


A  =  iw 


.  2 
V 


-I  2/  3 
-/A  x  , 


(1.5) 


The  path  of  the  outward  moving  electrons  is  obtained  by  integrating 
Eq.  (1.5),  namely, 


x  =  3 


Va  j^v 


3/2  ,  1/2  r-  3 

o  -  <v0  ~  v* 


(1.6) 


We  observe  that  the  time  between  emission  and  maximum  displacement  is 

(1.7) 


v  1/2 

•  (i)  ■ 


and  the  time  between  emission  and  return  is 


v  1/2 

t  ~t  =  2t  =  2  ( 
r  1  \  A  / 


The  total  number  of  electrons  outside  the  emitting  surface  is  then 

v  1/2 

=  n0(tr-tl)  =  2n0(^)  . 


(1.8) 


N 


Or 


N 


W) 


1/2 


ss 


(1.9) 


The  dipole  moment  is 
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Z(t)  = 

Substituting  from  Eqs 
extends  from  0  to  the 

Z(t)  = 


-e^dx  p(x)  x  . 

.  (1.2)  and  (1.5)  and  observing  that  the  integration 

3/2  1/2 

maximum  displacement  /  3A  ,  we  find 

2 

mv0 

-e  - g  . 

87re 


(1.  10) 
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Appendix  2 

SOLUTION  FOR  t  >  T 


For  a  constant  emission  rate  n^  and  using  the  no— screening  approxi¬ 
mation,  Eq.  (7)  gives  for  the  motion  of  an  individual  electron 

2 


4?re  -T 

x(t)  = - N  , 

m  0  ’ 


(2.1) 


where  Nq  =  n^r  and  r  is  the  time  interval  between  emission  and  return 
of  the  electron.  Thus  if  the  electron  is  emitted  at  t  =  t^  , 


x(t)  =  -  1§-  No{t_to)2  +v0(t"V- 


[2.2) 


From  this  result, 


mv0  N0 


2,re  No  “o 


Thus 


j  ui*0 

N„  ■  V"o  ■ 


(2.  3) 


(2.4) 


Assume  now  that  Eq.  (2.4)  is  valid  for  t  <  T  and  that  for  t  >  T, 
n(t)  =  n  e-^^  Then  from  Eqs.  (2)  and  (10), 


N(t)  =  (l  ~  ^  j  —  n^t^t)  ”  T)>  for  T  -  t  <  tj(t)  <  T, 

n„  ,  -/3(t,(t)-T) 


o  / 

r  (e 


-0(t-T) 


)  ,  for  t^t)  >  T  , 


(2.5) 


RM-2849 

29 


where  t^(t),  the  time  of  emission  of  an  electron  which  returns  at  time  t, 
is  given  by 


t  =  t  (t)  + 


2ffe 


1 

N(t^(t))  * 


(2.  0) 


Eqs.  (2.5)  and  (2.6)  can  be  solved  approximately  in  two  limits  - — 

£{t  —  t^(t))  <<  1  and  P(t~t^(t))  >>  1,  where  t  and  t^(t)  >>  T.  In  the 
case  that  (3(t  —  t^(t))  <<  1,  t  >>  T  and  t^(t)  >>  T,  Eq.  (2.5)  gives 

N(t)  %  n0eH3(t_T)  (t  -  t1(t»  .  (2.7) 


Substitution  of  Eq.  (2.  7)  into  Eq.  (2.6)  leads  to 


N(t) 


n# 


mv0no 


2jre 


(2.8) 


In  the  case  that  /3(t  —  t ^ (t>)  >>  1,  t  >>  T  and  t^(t)  >>  T,  Eq.  (2.5)  gives 


n  H3(t  ,(t)-T) 

N(t)  %  j-  e 


(2.9) 


Substitution  of  Eq.  (2.9)  into  Eq.  (2.  6)  gives 


M(t)  £ 


2^e2 


(2.  10) 


mV 

2- 

2 ire  n,. 


0(t  -  T)  +  In 
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Appendix  3 

INCOHERENT  RADIATION 


We  will  show  here  only  that  the  incoherent  part  of  the  radiated  field 
is  small  compared  to  the  coherent  part.  Eor  this  purpose  it  is  sufficient 
to  overestimate  the  intensity  of  the  incoherent  field.  For  an  overestimate, 
it  will  be  assumed  that  each  electron  radiates  independently  at  all  frequencies. 
The  no— screening  approximation  will  be  used  and  it  will  be  assumed  that  the 
electrons  are  in  a  steady  state  corresponding  to  a  constant  y— ray  flux  over 

Q 

a  time  1/a  &  10  sec.  In  this  case  the  total  power  per  unit  area  radiated 
to  the  distance  R  is 


87 r  2  e  -.2  1 

P*T  a  No  3  x  ~~2 

c  4wR 


(3.  1) 


where  x  is  the  acceleration  of  each  electron  as  given  by  Eq.  (2.  1)  and 


Nq  is  the  number  of  electrons  outside  the  emitting  surface  as  given  by 


Eq.  (2.4).  Then 


3  6  2  ,mv 

128tt  e  a 


2  3 
m  c 


(z*  %) 


3/2 


2,ire 


4ttR 


(3.  2) 


where,  using  Eq.  (1), 


n.S;  — 
0  2\ 


K  a 
e  7.  5  x  10 


2i 


4tra 


3 

Using  a  =  10  cm, 
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P  % 


7  x  10 


R 


v3/  2  ,2 

Y  ergs/ cm  sec. 


km 


(3.3) 


Thus 


<  E  >  8s 


5  x  10  2Y3/f4 


R 


volts  /meter  . 


km 


(3.4) 


Comparison  of  Eq.  (3.4)  with  Eq.  (28)  establishes  that  the  incoherent 
radiation  is  negligible. 
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Appendix  4 

NOISE  FLUCTUATIONS 

Cosmic  background  noise  is  generally  described  by  specifying  an 
equivalent  "noise  temperature"  for  the  frequency  band  of  interest.  This 
equivalent  temperature  is  related  to  the  average  cosmic  noise  power  per 
unit  area  per  cycle  per  second  at  wavelength  X  by 

kT  —2  2  T  o 

SN  =  4?r  ~ 2  =  ^  x  10  —watts/m  eps,  (4.1) 


where  T  is  in  degrees  Kelvin  and  X  in  meters. 

2 

If  a  receiver  with  an  antenna  area  A  m  and  band  width  cr  cps 
integrates  the  received  noise  power  for  a  period  of  t  seconds,  the  total 
received  energy  will  be,  on  the  average, 

_  —22  A 

e  =  S  err  A  =  1.  7  x  10  T  — r-  Cr  t-  joules.  (4.  2) 

X^ 

The  fluctuation  of  the  mean  noise  energy  is  given  by  the  fluctuation  of  the 
radiant  energy  contained  in  a  hohlraum  of  volume  Act  and  temperature  T. 

For  the  band  width  O',  thi3  fluctuation  in  received  energy  is 


Act 


„  ,2  4  hv/kT 
Birh  v  e 

3.  hv/kT  .. 
c  (e  -1) 


o_ 

2  • 


For  the  case  of  radio  waves,  hv  <<  kT,  and 


(4.  3) 
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=  Arc 


A- 4) 


Thus  the  dispersion  of  the  mean  received  noise  energy  is 


